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Entangled many body systems have recently at- 
tracted significant attention in various contexts. 
Among them, spin squeezed atoms and ions have 
raised interest in the field of precision measure- 
ments, as they allow to overcome quantum noise 
of uncorrelated particles [H, 3, 4, 5, 6]. Pre- 
cise quantum state engineering is also required 
as a resource for quantum computation, and spin 
squeezing can be used to create multi-partite en- 
tangled states j?!. Two-mode spin squeezed sys- 
tems [8] have been used for elementary quantum 
communication protocols [9, 10]. Until now spin 
squeezing has been always achieved via gener- 
ation of entanglement between diff"erent atoms 
of the ensemble [111, [H, |l3]. In this Letter, 
we demonstrate for the first time ensemble spin 
squeezing generated by engineering the quantum 
state of each individual atom. More specifically, 
we entangle the nuclear and electronic spins of 
10^^ Cesium atoms at room temperature. We ver- 
ify entanglement and ensemble spin squeezing by 
performing quantum tomography on the atomic 
state. 

Ensembles of neutral atoms in the form of cold gases 
as well as room-temperature vapor have emerged as im- 
portant systems for the generation of entan glem ent and 
for the storage of quantum information [ol. llOl. 
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161 . [171, [181, [19|, [20|. It has been shown that ensemble 
spin squeezing can be used to significantly improve the 
fidelity of some existing quantum memory protocols P|. 
In the majority of the work on ensemble entanglement 
alkah atoms used as memory units were treated as spin- 
1/2 systems, while their actual much higher total angular 
momenta were considered an unfortunate complication. 
Entanglement in those experiments could therefore be 
generated only between different atoms of the ensemble. 
In this Letter, we generate entanglement in an ensem- 
ble of spin-4 atoms by using their internal structure to 
squeeze the individual spins. Indeed, we can show that 
spin squeezing within the Cesium F — A manifold is a 
signature of entanglement between the electron and the 
nucleus. However, as discussed below, the ensemble spin 
squeezing does not necessarily follow from the single atom 
spin squeezing. Therefore, we verify the collective squeez- 
ing by measurement of the collective spin state via strong 
coupling to off-resonant light. In particular, we perform 
for the first time quantum tomography of a non-classical 
state of an atomic ensemble. 



We first consider a collection of ground state atoms, 
each prepared in the magnetic sublevel \F — A,mp = 4) 
with respect to the x quantization axis. This fully 
stretched state is a product state, where both the nu- 
clear and the electronic spin are in Eigenstates of their 
angular momentum projection operators with maximal 
Eigenvalues Ix — 7/2 and Sx — 1/2, respectively. In the 
following, we denote the total spin operator of the en- 
semble J. For large atom numbers N, the total spin 
becomes macroscopic with (J) = (J, 0, 0) in cartesian 
coordinates, and J = TiNF. The commutation rela- 
tion [Jy^Jz] — ihJx imposes an uncertainty relation, 
and the fully pumped ensemble is in a minimal uncer- 
tainty state with (AJy)'^ = (AJ^)^ = hJ/2. Such states 
are called coherent spin states (CSS) and set the stan- 
dard quantum limit (SQL) to spin projection measure- 
ments. For a squeezed spin state (SSS), one of the or- 
thogonal projections exhibits an uncertainty below the 
SQL. For useful spin squeezing, however, it is not suf- 
ficient to reduce the uncertainty below its initial value 
(X^ = 2(AJj_)2/?iJ < 1), but to compare it to a CSS 
with the same mean spin (C^ ~ 2(A J^)^/?i(| J|) < 1) 
[H im , or even more rigorously to the initial uncert ainty 
in spin angle (f^ = 2J{AJ±)^ /h{\3\)^ < 1) 0, 0- 
The last condition is the strongest and provides a suf- 
ficient condition for entanglement of elementary spin-1/2 
constituents 2j|. It has to be met to improve quantum 



metrology with a given number of atoms, e.g., the preci- 
sion of atomic clocks. 

Unconditional spin squeezing can be achieved by ap- 
plying a Hamiltonian that depends non-linearly on carte- 
sian spin components orthogonal to the mean spin 21 1. 



A Hamiltonian of the form Hg = aF^ acting on a spin 
F induces a process called one-axis twisting. Starting 
from a CSS aligned along the x-direction, the spin is ro- 
tated about the z-axis in proportion to its fluctuating 
z-component. The initially symmetric quasi-probability 
distribution (QPD) will be sheared and appear squeezed 
in the y,z plane in a slightly rotated basis Fy',Fzi. At 
the same time, the magnitude of the mean spin is reduced 
as the QPD bends around the spherical phase space for 
constant F. 

In our experiment we generate a squeezing Hamil- 
tonian using the tensor light shift induced by an off- 
resonant light field. In case of linear polarization along 
the z-direction with classical field amplitude , the in- 
teraction with an atomic spin F reduces to an atomic 
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FIG. 1: Pulse sequence and level scheme for atomic state 
preparation and quantum state reconstruction, (a) Optical 
pumping creates a macroscopic spin J in a CSS, aligned with 
the magnetic field B. (b) A Raman transition creates a coher- 
ent superposition of even magnetic sublevels, which shears the 
QPD in the rotating frame and decreases the mean spin, (c) 
The squeezed spin state of the ensemble is characterized via 
off-resonant Faraday interaction with a probe beam, propa- 
gating in the ^-direction. Quantum noise limited polarization 
measurements reveal the statistics of the atomic angular mo- 
mentum projections Jyi and J^/, which are precessing about 
the X-axis. Simultaneously, a weak measurement along the x- 
axis determines the magnitude of the mean atomic spin (Jx)- 
See methods section for details. 



Hamiltonian of the form: 
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While the scalar polarizability ao simply results in an 
overall Stark shift, its tensor part, proportional to a2, 
provides the basis for spin squeezing. 

In the x-basis, the z-polarized field can be decomposed 
into equally strong u+ and a_ components, driving Ra- 
man transitions between magnetic sublevels (Am = 2). 
In the presence of a magnetic field along the x-direction, 
the above description remains valid after transformation 
to a frame co- rotating with the Larmor frequency uj^. 
Consequently, the frequencies of the two Raman fields 
have to be shifted by ±u;l to meet the two-photon res- 
onance condition (see Fig. [IJb)). In this sense, spin 
squeezing corresponds to the creation of paired excita- 
tions of the spin at frequency and will reveal itself in 
the Fourier component of the rotating spin at frequency 

Naturally, this type of interaction is non-linear only 
in spin components of individual atoms. Therefore, the 
maximally possible degree of squeezing ff\ is limited by 
the total spin of a single atom, while the total spin of an 



ensemble can be arbitrarily high. With one-axis twist- 
ing, the squeezing is further limited by the deviation 
of the resulting QPD from a geodesic in the spherical 
phase space. This finding is equivalent to the differ- 
ential light shift imposed on magnetic sublevels in the 
x-basis, which inevitably accompanies the Raman cou- 
pling and detunes the coupling between different sub- 
levels. A somewhat higher degree of squeezing can be 
achieved with a two-axis count ert wist ing Hamiltonian of 
the form = a{F^ - F^)^. Taking the decreas- 
ing mean spin into account, squeezing parameters of 
~ 0.163, ~ 0.247, and ~ 0.327 can be reached 
with this Hamiltonian, using the maximum available spin 
of = 4 in ground state Cesium atoms. The required 
Hamiltonian could, e.g., be implemented with two laser 
fields of orthogonal polarizations (along rotating y'- and 
^'-directions) and opposite detunings from an atomic res- 
onance. Here, we balance the light shift instead with 
the second-order Zeeman shift Hz = PP^, such that 
all energy splittings between magnetic substates of the 
F — 4-manifold coincide. By adjusting the light intensity 
and detuning, a2\Ez\'^ = —8/3 can be chosen to obtain: 
Hs+Hz = -^^^ + ph^F{F + l)/2 + PiF^ - F2)/2, 
using F^ + F^ + F^ = h^F{F + 1) for Eigenstates of F^. 

Coherent superpositions of magnetic sublevels similar 
to those described above have been recently generated by 
universal quantum control of a hyperfine spin in an ultra- 
cold ensemble [25]. In that experiment the mean values 
of elements of a single atom density matrix were deter- 
mined. However, knowledge of the single atom density 
matrix is not sufficient to infer spin squeezing of the en- 
semble. Due to the collective preparation, the members 
of the ensemble will not in general be uncorrclated. Col- 
lective spin squeezing, as required for quantum limited 
metrology or quantum memory Q , can therefore only be 
determined via measurement of collective quantum fluc- 
tuations, ideally via full quantum tomography, as per- 
formed in the present paper. Only such measurements 
can reveal the deleterious effect of classical and quan- 
tum atom-atom correlations, which can be created in the 
preparation process and can easily ruin the spin squeez- 
ing of the ensemble. In particular, it is of major impor- 
tance to send the driving Raman field (Fig. [Hb)) along 
the direction of the mean spin (x-direction) . If, e.g., the 
orthogonal y direction is used, the quantum fluctuations 
of the driving field will couple via atomic 7r-transitions to 
the collective symmetric mode for Am^ = 1 coherences, 
and thus destroy collective spin squeezing under strong 
coupling conditions. In fact, this process forms the basis 
for quantum memory applications [lo| . 

The protocol for atomic quantum tomography used to 
reconstruct the collective quantum state is detailed else- 
where [27l . [2^. In brief, we send a strong, linearly 
y-polarized beam along the z-axis of the atomic ensem- 
ble, see Fig. [IJc). After the interaction, we perform po- 
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FIG. 2: Exemplified quantum state reconstruction (at « 3 dB 
squeezing) with ~ 0.8. (a) Scatter plot of 10* realizations. 
The scaled atomic covariance (red ellipse with v^c-radii) is 
determined by correcting the total covariance (blue) for light 
and back-action noise (green), (b) Maximum likelyhood es- 
timation of the density matrix p for the collective spin state, 
showing alternating excitation numbers m, n in the harmonic 
oscillator approximation. See methods section. 
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FIG. 3: (a) Results of the quantum state reconstruction for 
varied Raman pulse length. Red and blue diamonds show 
atomic squeezed (= x'^/2) and anti-squeezed variances, re- 
spectively. For useful squeezing, they have to be compared 
to the variance of a CSS with the same mean spin (black 
squares), (b) Resulting squeezing parameters (diamonds) 
and (crosses), representing renormalized data from the in- 
dicated region in (a). Theoretical predictions are shown as 
solid and dashed lines. 



larization homodyning in the 45°-basis, measuring the 
photon number difference in the two outputs, i.e. the 
Stokes operator Sy = ri+450 — n_45o . Due to the Faraday 
effect, the polarization of the input beam is rotated pro- 
portionally to the instantaneous atomic spin component 
Jz, producing a proportional homodyne output for small 
rotations. The atomic spin precesses, and information on 
the two orthogonal, rotating spin components Jyi and Jz' 
can be obtained by evaluating the cosine and sine com- 
ponents of the output signal at the Larmor precession 
frequency lol- The atomic ensemble is freshly prepared 
and the procedure repeated 10^ times to determine the 
quantum statistics including the mean and variance of 
the collective atomic spin state. 

The state tomography is illustrated in Fig. [2l Atomic 
variances are directly inferred from the sampled data. In 
addition, we performed maximum likelyhood estimations 



of the collective state in the 10 best-resolved dimensions 
of Hilbert space [2^. The results, i.e. the measured vari- 
ances, are shown in Fig.[3]together with predictions from 
a density matrix calculation. Our model includes atomic 
decay due to spontaneous emission and incorporates in- 
dependently measured values for depolarizing and de- 
phasing times in the dark (Ti = 80 ms, T2 — 20 ms). We 
kept the effective light power in the cell as a free model 
parameter, but especially for short times, the theoreti- 
cal results are stable against small variations in the con- 
trol parameters. The used light power slightly overcom- 
pensates the second-order Zeeman-shift, and our model 
suggests better performance at even higher powers as it 
identifies atomic decay as the dominant limitation. We 
consistently achieve squeezing of a canonical variable by 
— 0.47 ±0.1. Quantum memory [§] and metrology 
applications gain by = 0.54 ± 0.1. 

In conclusion, we demonstrated unconditional squeez- 
ing of a collective spin with a noise reduction of « —3dB 
by generation of entanglement within individual mem- 
bers of the ensemble. We have fully characterized the 
non-classical collective atomic state via quantum state 
tomography. The method developed in this paper is in- 
dependent and complementary to the ensemble entan- 
glement and squeezing approach based on quantum- non- 
demolition (QND) measurements [^, [l^, [lj|. Indeed it 
can, e.g., be directly applied, either separately or in com- 
bination with the QND method, to enhance the fidelity 
of the quantum memory [9(] ■ It may also find use in spec- 
troscopic applications, such as magnetometry beyond the 
SQL. 



METHODS 
State preparation 

In our experiment, we employ a room temperature 
ensemble of iV « 10^^ ^■^'^Cs atoms, contained in a 
parafhn coated, 22 x 22 x 22 mm'^ glass cell. The cell 
is shielded against field fluctuations and subjected to 
a homogeneous magnetic field along the x-direction of 
Bx ~ 0.9 G. By fine-tuning the magnetic field we adjust 
the atomic Larmor precession to a reference frame rotat- 
ing at ojl ^ 2ttx 322 kHz, which is defined by a computer 
controlled radio frequency (rf) synthesizer. 

The experimental pulse sequence consists of three 
stages and is depicted in Fig. [T] The two-step state 
preparation starts with an 8 ms long optical pump pulse 
of two resonant, circularly polarized beams, propagating 
along the a;-direction. See Fig. [2a). After optical pump- 
ing, we find 98% of the atoms in the \F = 4, nix = 4) 
state, by evaluating the magneto-optical resonance sig- 
nal (MORS) [2^. At the second stage, we coherently 
transfer atomic population between even magnetic sub- 
states of the F = 4-manifold by driving Raman transi- 
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tions on the Di-line for a variable duratfon T^, = — 6 ms 
with a single photon detuning of = —550 MHz from 
the Q'^Si/2,F = 4 ^ &'^Pi/2,F' = 4 transition. See 
Fig. mb). The two-photon difference vb. ~ 644.2 kHz 
is resonant with the nix = 4 — > nix = 2 transition, 
including corrections for second-order Zeeman and dif- 
ferential light shifts. Finally, the collective atomic state 
is analyzed via off-resonant Faraday interaction with a 
strong (P « 5 mW), linearly polarized, 2 ms long probe 
pulse. The probe laser frequency is close to the D2- 
line with a blue-detuning of Ap = 825 MHz from the 
F = 4 ^ F' = 5 transition. See Fig. ^c). 

Some experimental issues have to be considered for the 
generation of collective squeezed states. The state prepa- 
ration must be completed on a timescale short compared 
to decoherence mechanisms, e.g. caused by inhomoge- 
neous fields, atom-atom, and atom-wall collisions. There- 
fore, the strength of the squeezing interaction should be 
maximized. Although atomic motion leads to averaging, 
it is necessary, particularly on short time scales, to pro- 
vide rather homogeneous coupling strength over the cell 
volume with expanded, mode-matched Raman fields. A 
small angle between the two beams leads to a varying 
relative phase over the cell volume, which gives rise to 
a residual Doppler shift for moving atoms and results 
in incoherent squeezing. A similar effect is caused by 
the collective spin of the polarized ensemble. Due to 
the Faraday-effect, the sample is circularly birefringent 
and induces a phase shift over the cell length, A(/) < 12° 
for our parameters. A simple and robust experimental 
approach is to generate the necessary Raman fields by 
amplitude modulation of a single, linearly polarized light 
beam. But to avoid the additional sidebands and min- 
imize spontaneous emission, we mode-match two phase- 
stable, circularly polarized fields. For ground state alkali 
atoms, the tensor polarizability compared to the scatter- 
ing rate is limited by the excited state hyperfine splitting. 
Their ratio is optimal when tuned between the hyperfine 
transitions of the Di hne, (F = 4 ^ F' = 3, 4) at 894 nm 
for Cesium. At this detuning, the tensor polarizability 
is also largest and stationary and thus independent of 
atomic velocity. In addition, it has the correct sign for 
compensation of the second-order Zeeman shift. 



State reconstruction 

The commutation relation prevents us from measur- 
ing both spin components with arbitrary precision, and 
their precise values are masked by back-action of the 
original light noise. Their statistical moments, however, 
can be inferred from the moments of the light output. 
For large atom numbers, we can make a harmonic os- 
cillator approximation and define canonical operators as 

X = Jy' I \lfi{Jx), p = Jz' /y fi{Jx) for atomic variables 



and y ~ Sy/ y (Sx), q = 5'z/ \J (Sx) for light variables of 
a given mode with [x,p\ w [y, q] w i. The two simultane- 
ously accessible output variances are given by: 

(Ay^oUr = (Aj/D' + Y(A-i, Ap)2 + ^{Aq^r, (2) 

where the indices correspond to different, non-orthogonal 
temporal modes of the light field, oscillating as cosine (c) 
and sine (s) components at the Larmor frequency. Index 
refers to a rectangular envelope of duration T, and 
oc Jx T. 

We scale the time-integrated homodyne output by 
referencing it to the light vacuum noise, ((Ay)^) — 
((Ag)^) — 1/2. To avoid changes in the optical and 
electronic paths, this is measured by shifting the atomic 
Larmor frequency well outside the detection bandwidth 
using the magnetic field. The coupling strength is 
calibrated by measuring the resonant noise of the atomic 
ensemble in thermal equilibrium. For the thermal state, 
the back-action noise is zero and the output noise is 
given by the light noise and the atomic contributions 
from thermally populated states in the F = 4-manifold: 

principle, the second-order Zeeman-shift of the magnetic 
sub-levels has to be considered, but it is compensated 
by the differential light shift induced by the probe laser. 
The measurement of the mean spin Jx/N is calibrated 
with a fully pumped ensemble. 

The density matrix for the collective spin state is rep- 
resented in a basis given by excitation operators defined 
as b, — {x ± ip) /\/2 and m,n — b^b. See Fig. 
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